Directed Percolation and Generalized Friendly Walkers 



PL, 



John Cardy 1 ' 2 and Francesca Colaiori 1 
1 University of Oxford, Department of Physics - Theoretical Physics, 1 Keble Road, Oxford 0X1 3NP, U.K. 

2 All Souls College, Oxford. 

We show that the problem of directed percolation on an arbitrary lattice is equivalent to the 
problem of m directed random walkers with rather general attractive interactions, when suitably 
continued to m — 0. In 1+1 dimensions, this is dual to a model of interacting steps on a vicinal 
surface. A similar correspondence with interacting self-avoiding walks is constructed for isotropic 
percolation. 
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The problem of directed percolation (DP), first intro- 
duced by Broadbent and Hammersley |jj , continues to at- 
tract interest even though it has so far defied all attempts 
at an exact solution, even in two dimensions. Although 
the problem was originally formulated statically on a lat- 
tice with a preferred direction, when the latter is inter- 
preted as time the universal behavior close to the percola- 
tion threshold is also believed to describe the transition 
from a noiseless absorbing state to a noisy, active one, 
which occurs in a wide class of stochastic processes M. 
It also maps onto reggeon field theory, which describes 
high-energy diffraction scattering in particle physics 11 . 

Some time ago, Arrowsmith, Mason and Essam Q ar- 
gued that the pair connectedness probability G(r,r') for 
directed bond percolation on a two-dimensional diagonal 
square lattice can be related to the partition function for 
the weighted paths of m 'friendly' walkers which all begin 
at r and end at r', when suitably continued to m = 0. 
These are directed random walks which may share bonds 
of the lattice but do not cross each other (see Fig. [I]). 
In fact, Arrowsmith et al H represented these configu- 
rations in other ways: either as vicious walkers, which 
never intersect, by moving the friendly walkers each one 
lattice spacing apart horizontally; or as integer flows on 
the directed lattice, to be defined explicitly below. Ar- 
rowsmith and Essam || showed that G(r, r 1 ) is also re- 
lated to a partition function for a A-state chiral Potts 
model on the dual lattice, on setting A = 1, thus general- 
izing the well-known result of Fortuin and Kasteleyn |fj] 
for ordinary percolation. 

In a more recent paper, Tsuchiya and Katori M have 
considered instead the order parameter of the DP prob- 
lem, and have shown that in d=2 it is related to a certain 
partition function of the same A = 1 chiral Potts model, 
and also that, for arbitrary A, the latter is equivalent to 
a partition function for m — (A — l)/2 friendly walkers. 

It is the purpose of this Letter to describe a broad gen- 
eralization of these results. We demonstrate, in particu- 
lar, a direct connection between a general connectedness 
function of DP and a corresponding partition function for 
m friendly walkers when continued to m = 0. We show 
that this holds on an arbitrary directed lattice in any 



number of dimensions, and for all variant models of DP, 
whether bond, site or correlated. Moreover the weights 
for a given number of walkers passing along a given bond 
or through a given site may be chosen in a remarkably 
arbitrary fashion, still yielding the same result at m = 0. 

We now describe the correspondence between these two 
problems in detail. A directed lattice is composed of a set 
of points in R d with a privileged coordinate t, which we 
may think of as time. Pairs of these sites (fi,rj) are 
connected by fixed bonds, oriented in the direction of 
increasing t, to form a directed lattice. In the directed 
bond problem, each bond is open with a probability p 
and closed with a probability 1—p, and in the site prob- 
lem it is the sites which have this property. In princi- 
ple the probabilities p could be inhomogeneous, and we 
could also consider site-bond percolation and situations 
in which different bonds and sites are correlated. Our 
general result applies to all these cases, but for clarity 
we shall restrict the argument to independent homoge- 
neous directed bond percolation. The pair connectedness 
G(r, r') is the probability that the points r and r' (with 
t < if) are connected by a continuous path of bonds, 
always following the direction of increasing t. 

On the same lattice, let us define the corresponding in- 
teger flow problem. Assign a non-negative integer- valued 
current n(ri,rj) to each bond, in such a way that it al- 
ways flows in the direction of increasing t, and is con- 
served at the vertices. At the point r there is a source of 
strength m > 1, and at r' a sink of the same strength. 
There is no flow at times earlier than that of r or later 
than that of r' . Such a configuration may be thought of 
as representing the worldlines of m particles, or walkers, 
where more than one walker may share the same bond. 
The configurations are labeled by distinct allowed val- 
ues of the n(ri, rj), so that they are counted in the same 
way as are those of identical bosons. Alternatively, in 
1 + 1 dimensions, we may regard the walkers as distinct 
but with worldlines which are not allowed to cross. In 
the partition sum, each bond is counted with a weight 
p{n{ri,Tj)). In the simplest case we take p(0) = 1 and 
p(n) = p for n > 1 (although we shall show later that 
this may be generalized). Since p > p n for n > 1, there 



is an effective attraction between the walkers, leading to 
the description 'friendly'. The partition function is then 

Z(r;r';m) = ^ II P( n ( r *,rj)) 

allowed configs (ri,rj) 

This expression is a polynomial in to and so may be eval- 
uated at m = 0. The statement of the correspondence 
between DP and the integer flow problem for the case of 
the pair connectedness is then 

G(r;r') = Z(r;r';0). 

Note that since the weights p(n) behave non-uniformly 
as n — > 0, the continuation of Z(r; r'; to) to to = is not 
simply the result of taking zero walkers (which would be 
Z = 1): rather it is the non-trivial answer G. Similar re- 
sults hold for more generalized connectivities. For exam- 
ple, if we have points (r[, r 2 , ■ ■ . , r[) all at the same time 
t' > t, we may consider the probability G(r; r[, r' 2 , ■ ■ ■ , r[) 
that all these points, irrespective of any others, are con- 
nected to r. The corresponding integer flow problem has 
a source of strength m > I at r, and sinks of arbitrary 
(but non-zero) strength at each point r'. In this case 

G{r-r' i y 2l ...y i ) = {-l) l - 1 Z{r-r' i y 2l ...y i -m^Q) 

where the partition function is defined with the same 
weights as before. Since the order parameter for DP may 
be defined as the limit as t' — t — > oo of P(i? — t), the prob- 
ability that any site at time t' is connected to r, and this 
may be written using an inclusion-exclusion argument as 

P(t> t) = $>(r;r') - £ G^r'^r',) + ■■■ 



(where the sums over the r' are all restricted to the fixed 
time t'), we see that it is in fact given by the to = 
evaluation of the partition function for all configurations 
of to walkers which begin at r and end at time t' . This 
generalizes the result of Tsuchiya and Katori |7|] to an 
arbitrary lattice. Although this continuation to to = is 
reminiscent of the replica trick, it is in fact quite differ- 
ent. Moreover it is mathematically well-defined, since, as 
we argue below, Z is a finite sum of terms, each of which, 
with the simple weights given above, is a polynomial in 
in. 

We now give a summary of the proof, which is ele- 
mentary. The connectedness function G(r; r[, r' 2 , ■ ■ . , r[) 
is given jq] by the weighted sum of all graphs Q which 
have the property that each vertex may be connected 
backwards to r and forwards to at least one of the r'. 
(Alternatively, Q is a union of directed paths from r to 
one of the r'.) Each such graph is weighted by a fac- 
tor p for each bond and (—1) for each closed loop. A 
simple example is shown in Fig. 0. A given graph corre- 
sponds to summing over all configurations in which the 



bonds in Q are open, irrespective of all other bonds in 
the lattice. The factors of (—1) are needed to eliminate 
double-counting. It is useful to decompose vertices in Q 
with coordination number > 3 by inserting permanently 
open bonds into them in such a way that the only vertices 
are those in which two directed bonds merge to form one 
(2 — ► 1), and vice versa. This does not affect the con- 
nectedness properties. We may then associate the factors 
of (—1) with each 1 — ► 2 vertex in Q, as long as we incor- 
porate an overall factor (— 1) /_1 in G. With each graph 
Q we associate a restricted set of integer flows, called 
proper flows, such that n > 1 for each bond in Q, and 
n = on each bond not in Q. Those corresponding to 
the graphs in Fig. || are shown in Fig. |[ Note that the 
last graph corresponds to TO — 1 configurations of integer 
flows, which gives precisely the required factor of (—1) 
when we set to = 0. In general, summing over all al- 
lowed integer flows will generate the sum over all allowed 
Q, with correct weights p: the non-trivial part is to show 
that we recover the correct factors of (—1) when we set 
to = 0. 

This follows from the following simple lemma: if 
A(n) is a polynomial in n, and we define the polyno- 
mial B(m) = Y,n=i A ( n ), then B{0) = -A{0). We 
give a proof which shows that the result may be gen- 
eralized to other functions: write A(n) as a Laplace 
transform A(n) = f c (ds/2m)e ns A(s). Then B(m) = 
J c (ds/2Tri)((e s - e ms )/{\ - e s ))A(s), so that 5(0) = 
— J c (ds/2m)A(s) = —A(0). An immediate corollary 
is that if A(ni,n 2 , . . .) is a polynomial in several vari- 
ables, and B(m) = Y^=i ^( n : to — n, . . .), then B(0) — 
— A(0, 0, . ..). We use this to proceed by induction on 
the number of 1 — > 2 vertices in Q. Beginning with the 
vertex which occurs at the earliest time, the contribu- 
tion to Z from the proper flow on Q, when evaluated at 
to = 0, is, apart from a factor (—1), equal to that for 
another graph Q 1 which will have one fewer 1 — > 2 ver- 
tex. However, Q' differs from the previously allowed set 
of graphs Q in that it may have more than one vertex 
at which current may flow into the graph. For this rea- 
son we extend the definition of the allowed set of graphs 
to include those in which every vertex is connected to 
at least one 'input' point (ri,r 2 ,...) and at least one 
'output' point (r[ 7 r' 2 , . . .). In the corresponding integer 
flow problem, currents (toi, m 2 , . . .) flow in at the inputs, 
whereas the only restriction on the outputs is that non- 
zero current should flow out. The partition function is 
then the weighted sum over all such allowed integer flows. 
Induction on the number of 1 — ► 2 vertices then shows 
that this partition function, evaluated at rrij = 0, gives 
the corresponding DP graph correctly weighted. (The in- 
duction starts from graphs with no 1 — > 2 vertices which 
involve no summations and for which the result is trivial.) 

Since our main result relies only on the lemma it fol- 
lows also for rather general weights p(n). The only re- 



quirement is that p{n) grow no faster than an exponential 
at large n, and that, when continued to n = 0, it give the 
value p =/= 1. In this case, Z will no longer be a polynomial 
in to, but, since by the inductive argument above it is 
given by a sum of convolutions ofp(n), its continuation to 
?7i = will be well-defined through its Laplace transform 
representation. For example, we could take p(n) = p 1 ^ 71 
for ?i > 1. This raises the possibility of choosing some 
suitable set of weights for which the integer flow prob- 
lem, at least in 1 + 1 dimensions, is integrable, for ex- 
ample by Bethe ansatz methods. Unfortunately our 
results in this direction are, so far, negative. In the 
case of bond percolation on a diagonal square lattice let 
Z(xi, X2, ■ • ■ , x m ; t) be the partition function under the 
constraint that the walkers arrive at {xi, x<i, . . . , x m } at 
time t, the physical region being {x\ < X2 < • ■ ■ < x m }- 
Turning the master equation for Z in an eigenvalue prob- 
lem M and writing the eigenfunction ip m (xi, x%, . . . x m ) 

in the usual Bethe ansatz form, one gets for ^{xi, x<z) = 
Al2 ^(x 1 k 1 +x 2 k 2 ) + A2ie i(x 1 k 2 +x 2 k 1 ) the f n owing con di_ 

tion on the amplitudes: 

A g i(fei-fc 2 ) _ £ ( e i{ki+k 2 ) _|_ g-i(fei+fe 2 ) s \ 

-4l2 e -i{ki-k 2 ) _ e (gi(ki+k 2 ) _|_ e -i(ki+k 2 )\ 

(the same as that which appears in the XXZ spin chain 
p0[) where e = p(2)/p(l) 2 — 1. Requiring that the m- 
particle scattering should factorise into a product of these 
two-body S'-matrices places constraints on the weights 
p(n). In general these equations appear too difficult to 
solve, except in the weak interaction limit (e ~ 1/2), 
where we find 

71-1 

2" = 2q(n) + y"q(n- s)q(s)(l - Xs(n - a)) + 0{\ 2 ) 



where q(s) = p(s) / (p(l)) s , A = 2e— 1. This may be solved 
for successive q(n), but it is easy to see by applying the 
above lemma that, when continued to n = 0, it will al- 
ways yield the value 1, rather than p as required. We 
conclude that the to = continuation of this integrable 
case does not correspond to DP. It is nevertheless inter- 
esting that integrable models of such interacting walkers 
can be formulated. 

In 1+1 dimensions, our generalized friendly walker 
model maps naturally onto a model of a step of total 
height 777, on a vicinal surface, by assigning integer height 
variables h(R) to the sites R of the dual lattice, such 
that h = for x — > —00, h = m for x — > +00, and 
h increases by unity every time the path of a walker is 
crossed. The weights for neighboring dual sites R and R' 
are p(h(R') — h(R)). This is slightly different from, and 
simpler than, the chiral Potts model studied in Mm . 

A similar correspondence between percolation and in- 
teracting random walks is valid also for the isotropic case. 
The pair connectedness G{r, r') may be represented by a 



sum of graphs Q, just as in DP ||. Each graph con- 
sists of a union of oriented paths from r to r', As before, 
each bond is counted with weight p and each loop car- 
ries a factor (—1). Note that graphs which contain a 
closed loop of oriented bonds are excluded. Such contri- 
butions cannot occur in DP because of the time-ordering. 
The correspondence with integer flows or friendly walk- 
ers follows as before. The latter picture is particularly 
simple. ?Ti walkers begin at r and end at r' . When two 
or more walkers occupy the same bond, they must flow 
parallel to each other. Since they cannot form closed 
loops, they are self-avoiding. Moreover, walkers other 
than those which begin and end at r and r', which could 
also form closed loops, are not allowed. Each occupied 
bond has weight p(n) as before, and the separate con- 
figurations are counted using Bose statistics. G(r,r') is 
then given by the continuation to 777, = of the parti- 
tion function. We conclude that ordinary percolation is 
equivalent to the continuation to to = of a problem 
of 777, oriented self- avoiding walks, with infinite repulsive 
interactions between anti-parallel segments on the same 
bond, but attractive parallel interactions. In two dimen- 
sions, this is again dual to an interesting height model, 
in which neighboring heights satisfy \h(R') — h(R)\ < m, 
but local maxima or minima of h{R) are excluded. For 
example, the order parameter of percolation is given by 
the continuation to to = of the partition function for a 
screw dislocation of strength 777 in this model. 

To summarize, we have shown that the DP problem 
is simply related to the integer flow problem, or equiva- 
lently that of to bosonic 'friendly' walkers, when suitably 
continued to 777, = 0. This holds on an arbitrary directed 
lattice in any number of dimensions, and with rather gen- 
eral weights. It is to be hoped that this correspondence 
might provide a new avenue of attack on the unsolved 
problem of directed percolation. 
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FIG. 1. Typical configuration of m — 4 friendly walkers on the diagonal square lattice. 



FIG. 2. Allowed graphs Q corresponding to G((0, 0); (0, 2)). The first two are counted with weight p . The last is necessary 
to avoid double-counting, and conies with weight — p . 
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FIG. 3. Sets of configurations of m friendly walkers corresponding to each of the graphs in Fig. 2. The last corresponds to 
m — 1 configurations. 



